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Abstract 

We analyze the interrelation between the coexistence curve of the Lennard-Jones fluid and the 
Ising model in two and three dimensions within the global isomorphism approach proposed earlier 
[V. L. Kulinskii, J. Phys. Chem. B 114 2852 (2010)]. In case of two dimensions we use the exact 
Onsager result to construct the binodal of the corresponding Lennard-Jones fluid and compare it 
with the results of the simulations. In the three dimensional case we use available numerical results 
for the Ising model for the corresponding mapping. The possibility to observe the singularity of 
the binodal diameter is discussed. 
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The connection between continuous and discrete systems in Statistical Physics is the 
source of new methods and fruitful applications. The lattice models are usually considered 
as the caricatural pictures of the real ones. Due to their numerical and analytical tractability 
they serve as the main sources of the rigorous results. Famous Onsager's solution of the two- 
dimensional Ising model is the well known example which gave new impact for the theory 
of Critical Phenomena. In the modern theory of the Critical Phenomena the ideology of the 
isomorphism classes of the critical behavior provides the description of the real systems using 
the results obtained for the model systems among which the lattice models play important 
role HJ. In particular the molecular liquids with short range interactions of the Lenard- Jones 
(LJ) type: 

m = m {-('-)\ {!)"), (!) 

belong to the isomorphism class of the Ising model or equivalently to the Lattice Gas (LG) 
model. The latter is determined by the Hamiltonian: 

H = - J) J UjUj - ^ Tij , (2) 

where = 0, 1 whether the site is empty or occupied correspondingly. The quantity J is 
the energy of the site-site interaction of the nearest sites i and j, \i is the chemical potential. 
The order parameter is the probability of occupation x — (tii) and serves as the analog of 
the density. The phase diagram of the LG is symmetrical with respect to the line Xq = 1/2 
and formally extends up to the region T — > 0, where the limiting states x = and x = 1 
exist only. 

In [[J the approach which extends the notion of the isomorphism between the LG and 
the LJ fluid from the critical region to the whole liquid- vapor part of the phase diagram was 
proposed. It is based on the fact that the line x = 1 of the LG can be thought of as the 
analog of the Zeno-line for liquid []3], |J. In construction of the mapping the results of works 
of Aphelbaum and Vorob'ev |5|, |(| on the Zeno-line was extensively used. 

The objective of this paper is to use the global isomorphism transformation to the con- 
struction of the binodal of the Lennard- Jones fluids. We relax the condition of the constancy 
for the compressibility factor used in above mentioned works of Aphelbaum and Vorob'ev 
and interchange less restrictive. We will demonstrate that in such case we are able to get 
rather good estimates for the locus of the critical points and map the binodals of the Ising 
model to that of the Lennard- Jones (LJ) fluid. 
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The following mapping representing the global isomorphism between the thermodynamic 
states of the LG and the LJ fluid was constructed: 



with 



T 



T z -T c 

Here n and T are the density and the temperature of the fluid, t is the temperature variable 
of the LG normalized to the critical temperature so that t c — 1 we also use the standard 
dimensionless values for T and n of the LJ fluid J7|. Tz and rib are the parameters of the 
linear element: 

- + ^ = l. (4) 
The coordinates of the CP for the fluid are: 

n c = — ; — - — - , T c = Tz . (5) 

2(1 + z) 1 + z K J 

In contrast to the approach of || the parameters Tz and rib are chosen from the condition 
of consistency with the van der Waals (vdW) approximation for the given equation of state 
(EoS). In particular, Tz equals to the Boyle temperature in the vdW approximation. The 
parameter rib is chosen from the condition of linear change 

d Z(n(T),T)-l 



dT n(T) 

of the compressibility factor Z: 



0, (6) 



Z = — = l- — r 3 —^g 2 (r;n,T)dr, 7 
nl 61 J or 

along the linear element (§]). Using the virial expansion it is easy to get the condition: 

nb - Tz B~(Tz-)' (8) 

which formally coincides with the standard one (see e.g. H). The difference is that we relax 
the condition of the Constance Z along the linear element (f|) since it is valid only for the 
vdW EoS. As was shown in || such definitions along with the condition 

l/z = l + s/d, (9) 
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FIG. 1. Fitting the binodal of two dimensional LJ fluid (pi) with z = 1/3, Tz = 2.03 and = 0.971 



to the results of the simulations [1C, 11]. Also the corresponding locus of the CP is shown. 



following from the scaling properties of the attractive potential $ 



attr 



-(s+d) 



give the 



estimates for the CP loci of 2D and 3D LJ fluids which are in good agreement with the 
results of the numerical simulations. 

To prove that such closeness is not a coincidence we check another consequence of (^]). 
The latter provides the mapping between the binodals of the LG (Ising model) and LJ fluid. 
In accordance with ([|) z = 1/3 for 2D LJ fluid and z = 1/2 for 3D case respectively. For 
2D case the known exact solution for 2D Ising model is: 



x = 1/2 ±f(ty/\ /(*) = !-- 



1 



sinh 4 (2 J ft 



(10) 



Substituting it into (|3]) the parametric representation for the binodal of LJ fluid is obtained 
and can be compared with the known numerical simulation results (see Fig. |l|). Note that 
the values of the parameters Tz and obtained by fitting are very close to that following 
from the consideration above Tz ~ Tg dW = 2 and nj, ~ 0.91 (see also ||). The locus of the 
CP obtained is in good correspondence with the results of 



The same procedure can be applied to the 3D case. Here we use the results of |12| for 
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FIG. 2. Binodal for 3D LJ fluid (blue line) obtained from (fid]) the mapping after transformation 



with z = 1/2, T z = T 



(vdW) 
B 



4. The grey dashed line shows the range of applicability of (|iT|). 



The red line is for the Guggenheim binodal (|12j). Circles are for the data [14, 15] (different circles 
at a given temperature in the gas branch are for different total number of particles in simulations 



of 0) 



the 3D Ising model. There data for the spontaneous magnetization M(t) = 2x(t) — 1 were 
described by the expression: 

1 



x = - (1 ±t p (6 + aiT + a Al T Al )) , T=l-t 



(11) 



where (3 = 0.326941, A x = 0.50842, b = 1.6919, a x = -0.42572366 , a Al = -0.34357731. 
This expression represents the data in the region 0.0005 < r < 0.26. So we use (0) as 
the "true" binodal of the 3D LG. From the other side for the real molecular fluids the 
Guggenheim's |L3| empirical expressions: 



n i,9 _ 1 , 3 M T/T ^ JL 7 fl T 1 IT 1 W3 



1 + - (1-T/T C )±- (1-T/T C ) J 



(12) 



are well known. 

In Fig. |] the result of transformation of the Ising binodal (pT|) is shown. The mapping of 



the Ising model quite satisfactory correlates with the data near the CP (r < 0.1), while the 
Guggenheim curve better reproduce the data in outer region far from the CP. 

Note that both (|Tl| ) and fll2|) do not account for the diameter singularity of the density 
TB|. The possibility to apply the transformation (13) directly to the results of simulations of 



LJ fluids is facilitated by the fact that it neglects the singular fluctuational corrections to 
the regular behavior of the diameter. According to the theoretical approaches ]TT|-|TP|] the 
diameter of the density of the molecular liquid has the following structure: 

n d = -l = D w \ T \* + A-a M 1 -" + D 1 |r| + . . . (13) 

To the best of our knowledge none of the leading singular terms t 13 or r 1_Q has never 
been discussed in computer simulations. According to [EDj the diameter anomaly is difficult 
to observe by molecular simulations. The data of computer simulations |15], pO|-F25| are 
claimed to be consistent with the law of rectilinear diameter and commonly, the locus of 
the CP is obtained by the extrapolation of the data from the two-phase region using the 
classical rectilinear diameter law. Then the value n c is greatly influenced by the form of 
such extrapolation in view of wide opening of the binodal. 

It seems evident that LJ fluid does not possess any particle-hole symmetry which lead to 
regular diameter and the density variable should display both anomalies in its diameter as it 



does in the case of real and model fluids [26]. The results obtained for real substances show 
that \Di_ a / D 2 p\ > 10 |26|, |27|. So at least |r| 1_a anomaly may be observed. This poses the 
problem as to the values of the critical amplitudes D 2 /3, Di_ a and their dependence on the 
microscopic interactions. This question needs further investigation. 

As a summary we have demonstrated the usefulness of the transformation ([3]) in order to 
get the binodal of LJ fluid. By this we also confirmed the adequacy of the construction of 
the linear element (|J) as the key point in construction of the isomorphism transformation 
(H). This gives the alternative to the commonly used Zeno-line based considerations [§], [|. 
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